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ABSTRACT 

We have observed an area of as 27 deg 2 to an rms noise level of ^ 0.2 mJy at 15.7 GHz, 
using the Arcminute Microkelvin Imager Large Array. These observations constitute the 
most sensitive radio-source survey of any extent (> 0.2 deg 2 ) above 1.4 GHz. This pa- 
per presents the techniques employed for observing, mapping and source extraction. We 
have used a systematic procedure for extracting information and producing source cata- 
logues, from maps with varying noise and Mv-coverage. We have performed simulations to 
test our mapping and source-extraction procedures, and developed methods for identifying 
extended, overlapping and spurious sources in noisy images. In an accompanying paper, 
lAMI Consortiu m: Davie s et al.l d2010t) . the first results from the 10C survey, including the 
deep 15.7-GHz source count, are presented. 



1 INTRODUCTION 



The 9C survey JWaldram et al.l2003L[2oToh mapped 29 deg 2 of sky 
to 5.5 mJy completeness at 15 GHz, in addition to several larger 
and shallower areas. The Ryle Telescope (RT), which carried out 
the survey, has subsequently been reconfigured and re-equipped 
to form the Large Array (LA) of the Arcminute Microkelvin Im- 
ager (AMI: I AMI Consortium: Zwart etal] [2008). As part of this 
metamorphosis, three of the RT's eight antennas were moved to 
the north of the (almost) east-west line on which they originally 
stood, providing the telescope with north-south baselines. In addi- 
tion, new front-end receivers and back-end electronics, including a 
new correlator, were installed. The result is a telescope with much 
larger bandwidth and improved flux-density sensitivity, allowing 
us to extend our 15-GHz-band survey work to much deeper flux- 
density levels. 



* We request that any reference to this paper cites 'AMI Consortium: 

Franzen et al. 2010' 

t Email: t.franzen@mrao.cam.ac.uk 

$ Email: m.davies@mrao.cam.ac.uk 



The 9C survey was conceived to provide information regard- 
ing the foreground radio sources that contaminated the Cosmic Mi- 
cro wave Background ra diation observations of the Very Small Ar- 
ray dWatson et al . 2003). Similarly, the IOC source survey has been 
designed to complement the other AMI science programmes, which 
also require knowledge of contaminating radio sources. 

This paper is focussed on the techniques employed for ob- 
serving, mapping a nd source extraction in the IOC surv ey. In an ac- 
companying paper JAMI Consortium: Davies et alj20l61 hereafter 
Paper II) the first results from the 10C survey, including the deep 
15.7-GHz source count, are presented. 



2 THE ARCMINUTE MICROKELVIN IMAGER LARGE 
ARRAY 

The AMI consists of two separate telescopes - the LA and Small 
Array (SA). The IOC-survey data were collected using only the LA. 
Before proceeding further, the LA's vital statistics are summarised. 

The LA is an interferometer comprising eight 13-metre- 
diameter, equatorially-mounted dishes, with a range of baselines of 
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18-110 m. It operates at frequencies between 13.9 and 18.2 GHz 
with the passband divided into six channels of 0.72-GHz band- 
width. It has a primary beam at « 15.7 GHz of « 5.5 arcmin 
full width at half-maximum (FWHM) and a typical resolution of 
~ 30 arcsec (this varies depending on the precise Mv-coverage of 
any observation). The telescope measures a single, linear polarisa- 
tion (Stokes I +Q) and has a flux-density sensitivity of « 3 mjy for 
an integration time of one second. 

Accurate pointing is important for high-frequency observa- 
tions because of the relatively small telescope primary beams for 
such observations. The LA pointing is calibrated using five-point 
observations of bright point sources, which are carried out on a 
monthly basis. The data from five-point observations are used to 
construct a pointing model which provides the HA and Dec. offset 
of each antenna at any point in the sky. After applying these empir- 
ical corrections to the data, there is a residual rms antenna pointing 
error of « 30 arcsec. In practice, this is small enough compared 
with the primary beam size to have no significant effect on the ac- 
curacy of source flux densities measured from the raster maps. 



3 FLUX-DENSITY CALIBRATION 

It is standard practice, whilst observing using the LA, to visit a 
bright (> 200 mJy), unresolved source close (within a few degrees, 
if possible) to the target of interest for one in every 1 1 minutes. The 
majority of these sources are fou nd in the catalogue of the Jodrell- 
VLA Astrometic Survey (JV AS; |Patnaik et al.lll992l ; lBrowne et al.1 
1 19981 : 1 Wilkinson et alj|l998l) and are used for phase and amplitude 
calibration of the data. 

The LA is not yet sufficiently stable to rely on daily LA mea- 
surements of the primary calibrator sources for flux-density cali- 
bration, and so the secondary calibrators are used for flux-density 
calibration. The flux densities, in each channel, assumed for these 
secondary calibrator sources are obtained from SA measurements. 
The closest (in time) SA observation of the source is selected and 
the flux densities measured in each channel are adopted as the as- 
sumed values for calibrating the LA observation. As many of the 
secondary calibrator sources have varying flux densities, the SA 
observations are made within 10 days of the relevant LA observa- 
tion. In turn, the SA observations are calibrated by using SA obser- 
vations of the AMI primary calibration sources - 3C286 and 3C48 
- each of which are usually observed daily. 

The flux densities assumed for 3C286 and 3C48, in each of 
the AMI frequency channels, are shown in Table [T] The flux den- 
sities for 3C286 were converted from total-intensity measurements 
provided by R. Perley (private communication ) and made by the 
Very Large Array, and are consistent with the iRudv et al.l fl987h 
model of Mars transferred on to absolute scale, using results from 
the Wilkinson Microwave Anisotropy Probe. Values for 3C48 were 
obtained by measuring the ratios of the flux densities of 3C48 to 
3C286 using SA observations carried out between 2008 January 
and 2010 May. 

Finally, in order to account for phase errors, a correction factor 
of 1.082 is applied to raster flux densities. This correction factor 
was determined using self-calibrated maps of as 50 of the brightest 
sources in the survey fields - the data were collected using pointed 
observations, which were carried out to check the raster maps' flux- 
density scale and are described in Paper II. 



Table 1. Assumed flux densities for sources used for primary flux-density 
calibration (channels one and two are not used routinely, because of inter- 
ference from satellites). 
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Figure 1. Noise map for one of the 10C fields. Lighter shades indicate lower 
noise areas. 

4 OBSERVATIONS AND RASTERING TECHNIQUES 

The 10C survey covers a total area of « 27 deg 2 divided into 
10 fields. The area within the FWHM of the LA primary beam 
(~ 0.005 deg 2 ) is a small fraction of the total survey area. As a 
result, a rastering technique, similar to that used for the 9C survey, 
has been used. 

Each survey field is observed using a set of telescope point- 
ings which lie on a 2-D hexagonally-gridded lattice, projected on 
to the plane of the sky. The fields are observed with the lattice rows 
running along lines of constant declination at date. Since the tele- 
scope is equatorially mounted, this is the natural way to observe in 
order to minimise time lost to slewing between pointings. 

Maps are created for each of the pointings and combined 
as described below. No attempt has been made to recover spatial 
scales larger than the primary beam by jointly deconvolving data 
from separate pointings and, in this respect, the mosaicing tech- 
nique employed h ere is different from that of lCornwelil Jl988h and 
ISault etalJ i fl996l) . 

The spacing between pointing centres was chosen as 4 arcmin, 
as this was found (initially from simulations, but also in practice) to 
provide an acceptable compromise between the desire to minimise 
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variations in sensitivity across the final map (achieved by increas- 
ing the number of pointing directions) and to minimise observing 
time lost to telescope slewing. 

The spacing between pointing centres chosen for the IOC sur- 
vey (4.0 arcmin) is smaller than for the 9C survey (5.0 arcmin) 
in consequence of the LA's higher maximum observing frequency 
compared to the RT. Had the same spacing been used the varia- 
tion in sensitivity across the raster maps at the high-frequency-end 
of the AMI observing band would have been unacceptably greater 
than that on the RT raster maps (owing to the scaling of the primary 
beam with frequency). 

The time spent observing each field has been apportioned so as 
to create two distinct levels: the 'central', with lower noise, and the 
'outer', with a noise level about twice that of the central level. This 
scheme dovetails well with the SA observing programme. Survey- 
ing the outer areas to higher flux-density levels is also advantageous 
from the point of view of improving the source-count statistics at 
the brighter end of the survey. Fig. Q] shows the noise map for one 
of the 10C survey fields. The outer and central areas are readily 
apparent. 

Each of the fields was divided into smaller 'sub-fields', which 
come in two varieties: 'extended' which, taken together, cover the 
entirety of the field and 'core' which lie on the central area. Since 
the fields have different sizes, the number of pointings within each 
sub-field and the number of sub-fields vary between fields - the 
core sub-fields range between 54 and 72, and the extended between 
1 10 and 200, pointings. 

The sub-fields were each observed several times, typically for 
as 10 h at a time, until the desired noise levels were reached. The 
telescope dwell time per pointing was 30 s for the extended and 
60 s for the core sub-fields. The dwell times were chosen so that 
in an observing run of typical length, each pointing was visited 
several times, serving to improve the uv-coverage. Fig.[2]shows the 
uv-coverage for a typical pointing in a final raster map. Owing to 
flagging of the data, for example due to poor weather conditions 
during observing, the noise was sometimes found to be higher than 
desired over small areas of the final map; these areas were targeted 
with additional, small rasters or with individual pointings to ensure 
more uniform coverage over the field. 



5 MAPPING 

Raw data files from individual observing runs had various checks 
and calibrations applied using in-house, data-reduction software, 
reduce. A fuller explanation of the steps carried out in reduce, in- 
cluding the weig hting of data depending on the weather conditions, 
can be found in lAMI Consortium: Davies et alj d2009b . The data 
were then output as multisource wv-ftts files, which were, in turn, 
combined into a single multisource mv-fits file per field, again using 
in-house software. 

Depending on the precise nature of the flagging applied, it is 
possible that data belonging to individual telescope pointings are 
constituted from different frequency channels in slightly varying 
proportions. As a consequence, when the data were combined into 
a single wv-fits file, the data weights were adjusted, on a pointing 
by pointing basis, so that the continuum maps would always have 
the same centre frequency: 15.7 GHz. This causes a small loss of 
sensitivity. 
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Figure 2. The uv-coverage for a typical pointing in the core area of the 
final raster; data from all individual observations of this pointing have been 
combined. 

The mv-fits data were exported into the aip|3 package and in- 
dividual, continuum maps created for each telescope pointing using 
the imagr task. Each component map is 512 by 512 pixels, where 
the pixels are 5x5 arcsec 2 in size. Natural weighting was used to 
maximise signal-to-noise. 

Each component map was CLEANed using an elliptical Gaus- 
sian fitted to the central region of the dirty beam as the restoring 
beam. As a result, the restoring beam for each component map is 
slightly different - the implications of this for source extraction are 
discussed in Section [63] CLEANing was stopped after the first neg- 
ative clean component was reached, as real negative features are 
not expected to be present in the 10C maps. Negative clean com- 
ponents are in some cases necessary to best recover the true sky 
brightness distribution. However, we have chosen to adopt a con- 
servative approach in order to avoid the so-called 'clean bias' (aips 
cookbook; www.aoc.nrao.edu/aips/cook.html), where real source 
flux is underestimated as a result of attributing clean components 
to noise fluctuations. The noise on each component map was then 
estimated using the mean task, which fits a Gaussian, centred on 
zero, to the distribution of pixel values. 

Following a method similar to that of IWaldram et al.1 d2003h . 
the component maps were combined to form a single 'raster' map 
for each field, using the flatn task. Only data above the 0. 1 power- 
point of the primary beam from each component map were used. 
The mapping was also attempted using a cut-off of 0.3, since the 
uncertainties in the primary beam increase with distance from its 
centre. However, this was found to produce significant discontinu- 
ities in the map, which in turn caused spurious source detections 
(see Section [8j. Using a cut-off of 0.1 is found to make negligi- 
ble difference to the source flux-densities and to the map noise and 
avoids the problem of sharp discontinuities in the map. For any 
point on the final raster map, data from between one and seven 
component maps contribute. 

The map value M r at any point on the raster map is derived 
from the individual map, m,-, primary beam, p,, and noise, cr,-, val- 

1 astronomical image processing system - www.aips.nrao.edu/ 
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Figure 3. Detail of noise map for one of the IOC fields. Lighter shades 
indicate lower noise areas. The scale varies between 96 and 116 //Jy. 



ues of the relevant overlapping constituent maps, at that point. The 
maps are added, correcting for the primary beam, j- , and weighted 
according to the noise, (— ) 2 , such that 



Z mPi y IPi_ 



where j 6 (1,2, ...,7j. 



A noise map, which provides an estimate of the noise at the 
same point on the raster map (M„) is also computed, such that 



M n 



-1/2 



, where j 6 {1, 2, 7). 



Pixel values on both sky and noise maps are in Janskys per beam. 
The noise maps are found to provide an accurate representation of 
the noise. In particular, they are found to provide a better estimate 
of the true thermal noise close to sources, than would be found by 
considering the pixel values in the vicinity of sources on the real 
maps. Noise estimates computed using the latter method tend to be 
biased high. 

Given the 4.0-arcmin spacing of the pointing centres and the 
0.1 cutoff of the primary beam, if the noise level, cr c , was the same 
on each component map then the noise in the main area of the raster 
map would vary between » 0.86cr c and » 0.90cr c . Fig.[3]is a small 
section of the noise map for one of the survey fields and shows the 
structure of the noise. There is necessarily a sharp increase in the 
noise at the edges of the map, as seen in Fig.Q] 



6 SOURCE EXTRACTION 

Sources are detected in each survey field through analysis of the 
relevant raster map. For the 9C survey, candidate sources were fol- 
lowed up with pointed observations to allow identification of gen- 
uine sources. Observational time constraints meant that it was not 
possible to follow up each 10C source with a dedicated pointed 
observation. However, the IOC-survey noise maps provide much 
better estimates of the noise variation across the survey maps than 
were available as part of the 9C survey. Therefore, the noise maps 
were used to assess the reliability of each detection. Nevertheless, 
as described in Paper II, pointed observations were carried out to- 
wards as 50 of the brightest sources detected as part of the survey, 
providing a vital check of the raster maps' flux-density scale. 

In this section the steps in building a source catalogue are de- 
scribed. The methods employed are able to deal with noise levels 



and synthesised beams which vary across the maps. The source ex- 
traction is carried out using a combination of in-house software and 
aips. At the end of each subsection the parameters entered into the 
catalogue, at each of these steps, are provided. Source positions are 
quoted using equatorial coordinates (J2000). 



6.1 Source identification and peak measurement 

Source finding is carried out for each survey field making use of 
the noise map, which allows sources to be identified on the basis 
of their signal-to-noise ratios (SNRs). If one wishes to search for 
sources above ycr„, where cr„ is the value of the noise map at the 
pixel position and y is a constant usually chosen to be 5.0, an ini- 
tial search for pixels greater than 0.6ycr„ is conducted. The value 
of 0.6 was found to be sufficiently low so as to identify all peaks 
that have values greater than ycr„ after pixel interpolation. Pixels 
selected must have values greater than or equal to the surround- 
ing eight pixels. Pixels at the map boundary are not considered as 
potential sources. 

Since the synthesised beam is fully sampled (~ 6 pixels per 
FWHM of the beam), for each candidate source a peak value, cor- 
responding to a position interpolated between the grid points, is 
calculated. This is done by calculating the local map values on a 
successively finer grid (up to 128 times finer), by repeated convo- 
lution with a Gaussian-graded sine function of the form 



fix) 



^exp(-0.275x 2 ) 



if \x\ < 3 
otherwise, 



(1) 



where x has units of pixels dReedfl99"ol) . Interpolated peak values 
that are greater than or equal to y<x„ are stored. 

At this stage, the entries in the catalogue are S p k (the inter- 
polated peak flux density), and a p i and S p k (the interpolated peak 
position). The errors, SSp^, on the source peak- flux-densities are 
also provided. They are assumed to consist of a 5 per cent cali- 
bration error added in quadrature with the thermal noise estimate, 
because these are uncorrelated. Thus, the total error is SS p t = 
-\/o~l + (0.055 p k) 2 . The source names contained in the catalogue are 
constructed from their respective coordinates. 



6.2 Identifying overlapping/complex sources 

Following a method similar to that of Waldram & Rilevl Jl993l) . in- 
house software is used to identify overlapping sources. The integra- 
tion area of a source is taken to consist of contiguous pixels down 
to a lowest contour level of 2.5tr n . We begin by measuring the inte- 
gration area of the source with the lowest value of tr n . The source is 
classified as 'single' if its integration area contains no other peaks. 
Where more than one peak lies within the integration area, they are 
all classified as 'overlapping' and the name of the brightest peak, as 
well as the number of peaks, contained within the area are inserted 
in the 'group' column. 

The above procedure is repeated for each of the sources that 
remain unclassified, in order of increasing tr n . Considering the 
sources in this order avoids sources being classified as single when 
they actually lie within the 2.5cr n contour of another source. 

Sources can be separated by several synthesised beam widths, 
yet still share the same integration area because of the presence 
of faint extended emission linking them. Overlapping sources are 
therefore a good indicator of complex structure as well as point 
sources close together on the sky. 
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6.3 Gaussian-fitting 

Next, the aips task, jmfit is used to fit a 2D elliptical Gaussian, 
by an iterative least-squares process, to each source. Angular sizes 
and integrated flux densities can be estimated from these Gaussian 
fits. For each source, the fitted Gaussian is compared to the point- 
source response. This enables the morphology of resolved sources 
to be characterised. This also provides a check for 'over-narrow' 
sources, which appear significantly smaller than the point-source 
response and are likely to be spurious. Such 'sources' can be caused 
by discontinuities in the ma p, for ins t ance a t the boundary of a con- 
stituent map. According to Con don! i 19971) they can also occur at 
the positions of noise bumps and between the peaks of overlapping 
or complex sources. 

The comparison is complicated by the fact that the restoring 
beam is slightly different for each component map, as explained in 
Section [5] As most points on the raster map result from a combi- 
nation of several constituent maps, the point-source response varies 
across the raster and cannot, in general, be characterised by a single 
2-D elliptical Gaussian. In principle, this problem could be avoided 
by regrading the constituent maps with a single beam, having di- 
mensions that encompass all the individual restoring beams. How- 
ever, since in practice the beam areas vary significantly - by up 
to a factor of 2.24 for the most extreme field - this approach has 
not been adopted, as this would lead to a substantial loss of reso- 
lution. Rather, the point-source response belonging to the pointing 
with the highest weight at the position of the source is used in the 
comparison, as this proves to be a good approximation for the great 
majority of sources. 

Due to the large number of sources detected in the IOC sur- 
vey, the Gaussian-fitting has been automated. Six parameters are 
fitted for each single source: the peak value, peak RA and Dec, 
major axis, minor axis and position angle of the major axis. The 
fitting area consists of pixels inside a square centred on the source 
(peak) with a half-width of six pixels. Given the pixel size, this is 
sufficiently large to encompass point sources (for point sources at 
low SNR the 2.5-cr,, integration area is not sufficiently large). An 
exception is made if part of the integration area lies outside the 
box, in which case the size of the box is increased such that it en- 
compasses the integration area, allowing correct fitting of extended 
sources. The fitted parameters are initialised as S pk , a p ^, 6 p ^, 6 ma j, 
b m i n and b g (b m .^, i ra ; n and b 6 are the parameters of the restoring 
beam belonging to the pointing with the highest weight at the posi- 
tion of the source). A maximum number of 25 iterations is allowed 
to reach convergence. Even at low SNR, in real maps very few 
sources did not converge within this number of iterations. Those 
sources that failed to converge were found to lie on the sidelobes 
of bright sources and were often overlapping - they were excluded 
from the catalogue. 

For overlapping sources of N components, N Gaussians are 
fitted simultaneously. The peak value, peak position, major axis, 
minor axis and position angle of the major axis of each component 
are allowed to vary. The fitting area consists of pixels within a box 
encompassing N squares with half-widths of six pixels, one centred 
on each source (peak). Again, an exception is made if part of the 
integration area lies outside the box, in which case the size of the 
box is increased such that it encompasses the integration area. The 
values of 5 p k, a p k, <5 p k> ^maj> ^min an d b g for each source are used to 
initialise the fit. A maximum number of 25N iterations is allowed 
to reach convergence. 

Providing that the Gaussian-fitting achieves convergence, an 
integrated flux density, 5i„, is obtained. This is given by the peak 



Table 2. Interpretation of the values of the major and minor axes after de- 
convolution. 
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height of the fitted Gaussian multiplied by the ratio of the area of 
the fitted Gaussian to that of the point-source response. For com- 
plex sources that are many times the size of the synthesised beam, 
a more accurate integrated flux density can be obtained simply 
by summing pixel values from the peak down to some SNR, nor- 
malising with respect to the beam. However, this method has not 
been implemented here, since there are no believable sources in the 
IOC survey for which the Gaussian-fitting does not converge. The 
Gaussian-fitting also provides useful additional information. 

In addition, a^ n an d 6j n , the fitted peak position is recorded. 
Following the method o f lCondonl [l997), who has derived approx- 
imate and semi-empirical expressions for Gaussian-fitting errors in 
images with correlated noise, jmfit also provides an estimate of the 
error, cr in , on S m , due to the thermal noise. The thermal noise, cr„, is 
used by the task in estimating cr m . Following Section RTTl the total 

error on S m is taken as SS m = -^<rf a + (0.055 ; n ) 2 . 

An estimate of the source size is obtained by deconvolving the 
point-source response from the fitted Gaussian (both represented by 
2-D elliptical Gaussians). The sky brightness distribution is repre- 
sented by the resulting 2-D elliptical Gaussian, the major and minor 
axes, e ma j and e mm , and position angle, e s , of which are measured. 
For a point source in the limit of an infinite SNR this Gaussian will 
converge to a delta function. However, in the presence of noise the 
major axis will tend to be biased high (have a value greater than 0) 
and the minor axis low. In practice, axes for which the fitted value 
is negative are set to zero, as negative physical source sizes are not 
plausible. Table [2] summarises the interpretation of the major and 
minor axes, after deconvolution. 

The results of the deconvolution can be used to classify each 
source as point-like or extended. The classification criteria are jus- 
tified fully in Section [7] Briefly, the source is classed as extended 
if e maj is greater than or equal to some critical value, e crit , given by 
Equation [4] 

The catalogue entries for the Gaussian-fitting are, therefore, 
a in , S in , 5 in , (55 in , e crit , e ma j, e m i n and e e . A flag, f, is also included 
to denote one of two source types: point-like (P) or extended (E). 
As an additional check, for each source, the integrated flux density 
is normalised using the weighted sum of the areas of the point- 
source responses which contribute at the position of the source. If 
the resultant integrated flux density lies outside the range (5,„ - 
SS /„) to (5,„ + 6S in ) then a '*' is inserted in the 'flag' column. This 
indicates that the results obtained from the Gaussian-fitting must 
be treated with caution, since the error from the approximation of 
the point-source response is significant. In practice, this has been 
found to be the case for less than 1 per cent of sources above 5tr 
detected in all survey fields. 
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Figure 4. Section of an actual map. The 1-cr map noise is as 50 /^Jy. Contour 
levels start at ±100 /iJy and increase at each level by a factor of V2. The 
FWHM of the restoring beam of one of the component maps is displayed in 
the bottom left corner of the image. Source positions (ffpk,<5pk) found with 
y = 5.0 are represented as crosses. 



Table 4. The parameters which appear in the source catalogues. The section 
in which the parameter is described is shown in brackets. 



Source 


10C source designation J2000 {6j] 


Group 


10C group designation J2000 j6~2l 


ffpk 


RA (peak), in h, min, s (J2000) {6T} 


#pk 


Dec. (peak), in °, ', " (J2000) (6l\ 


Spk 


Peak flux density, in mJy )6.U 


<SS p k 


Error on peak flux density, in mJy 16.11 


«m 


RA (fitted peak), in h, min, s {63} 


Sis, 


Dec. (fitted peak), in °, ', "<6~3l 


S in 


Integrated flux density, in mJy 16. 3t 


SS i„ 


Error on integrated flux density, in mJy 16.3) 


Peril 


Critical component size, in " 16.3) 


^maj 


Major axis after deconvolution, in " 16.31 


^min 


Minor axis after deconvolution, in "{63) 


e$ 


Position angle after deconvolution, in ° , 




measured from North through East 163t 


t 


Source type (P = point-like, E = extended) 1631 


Flag 


A star indicates that the approximation error for 




the point-source response is significant {6.3\ 



6.4 Example catalogue 

A contour plot of a small section of one of the survey fields, cho- 
sen since it provides a particularly interesting example, is shown in 
Fig. [4] A subset of the catalogue entries for the sources detected 
in this region are provided in Table[3] Table|4]lists the parameters 
which appear in the final source catalogues. 



7 CLASSIFYING SOURCES AS POINT-LIKE OR 
EXTENDED 

As noted in Section [531 criteria for classifying sources as point- 
like or extended have been developed. One approach would be to 
use the residuals from 3-parameter fits, where only the peak bright- 
ness, and peak RA and Dec. are allowed to vary, to identify ex- 
tended sources; for extended sources the fit is expected to be poor. 
However, even for point sources the quality of the fit will vary with 
SNR and quantifying this dependence is difficult. An empirical ap- 
proach might be possible. However, a method, based on theory, 
which makes use of the angular size after deconvolution has been 
found to provide a straightforward test applicable to a wide range 
of SNRs. 

Maps containing point sources and noise were simulated in 
order to investigate the distribution of e ma j that would be expected 
for point sources as a function of SNR. Simulated visibilities were 
generated for point sources using wv-coverage from real observa- 
tions of a survey field. After simulating the sources and noise in 
the nv-plane, the resultant multisource ;/v-fits file was mapped as 
described in Section [5] and source extraction was carried out as in 
Section|6] 

Point sources were simulated in the uv-plane using a real mul- 
tisource kv-fits file as a template. The i th visibility associated with 
the / h pointing, Vy, was set to 

N 

Vjj = ^Skexpcftjjk , (2) 

k=l 

where N is the total number of sources, 5^ is the primary-beam- 
corrected flux density of the k' h source, is a phase factor given 
by 

</>ijk = 2?r (uij sina jk cos S jk + v tj sin 6 jk ) , (3) 

Uij and Vjj are the Mv-coordinates of Vy, and ajk and 6j k are the RA 
and Dec. separations between the k lh source and / h pointing centre 
respectively. 

A source was only added to any individual pointing if it fell 
within 15 arcmin of the pointing centre. Beyond this distance its 
contribution was negligible owing to the fall-off in the primary 
beam. Gaussian noise was also simulated, such that the noise on 
each of the component maps was identical, by adding random num- 
bers drawn from a Gaussian distribution to the real and imaginary 
parts of each visibility. 

We note that point sources can alternatively be inserted into 
real raster maps by adding Vjj to real multisource uv data. Such 
simulations were carried out in order to assess the completeness of 
the 10C survey (see Paper II). 

Our source extraction techniques were applied to a map of 
3.2 deg 2 containing 780, 100-mJy, simulated point sources - one 
at each pointing centre. The rms noise on each component map 
was simulated to be 10 mJy per beam. Once the component maps 
had been added together, as described in Section [5] the resulting 
rms noise varied between ss 8.6 and 9 mJy per beam over the cen- 
tral area of the raster map. For illustration, a section of this map is 
shown in Fig. [5] 

Fig. |6] shows the distribution of e ma j for the simulated point 
sources with SNRs of approximately 10. It can be seen that it is 
very rare for a point source to have e ma j > 40 arcsec at this SNR. 
However, as the SNR decreases, the distribution was found by do- 
ing runs with different source flux densities to move further away 
from zero and become wider and, as a result, 40 arcsec would not 
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Table 3. Abbreviated catalogue entries for sources detected in a section of an actual map (see Fig. [4). The flux-density correction accounting for phase errors, 
discussed in Paper II, has been applied. 



Source Group S pk (mjy) S in (mJy) e crit (") e maj (") e rain (") e fl (°) t 



10CJ002209+3 12803 




0.37 ± 0.06 


0.37 ± 0.10 


48.9 


9.0 




P 


10CJ002212+312317 




1.17 ±0.09 


1.13 ± 0.12 


27.5 


8.6 




P 


10CJ002223+312509 


100 002224+312409(3) 


1.46 ±0.10 


1.36 ±0.12 


25.0 


5.3 




P 


10CJ002224+312409 


10CJ002224+3 12409(3) 


6.9 ± 0.4 


7.1 ± 0.4 


25.0 


9.1 




P 


10CJ002227+312520 


10CJ002224+3 12409(3) 


0.66 ± 0.07 


0.59 ±0.10 


36.5 


4.3 




P 


10CJ002238+312622 




0.88 ± 0.08 


0.84 ±0.11 


31.0 


7.1 




P 


10CJ002246+313113 




3.33 ±0.19 


5.2 ±0.3 


25.0 


42.9 


9.0 


145.3 E 




. CD 



RIGHT ASCENSION (J2000) 



Figure 5. Portion of a simulated raster map with 100-mJy point sources 
lying at each pointing centre and noise about 10 mJy. Contour levels start 
at ±20 mJy and increase at each level by a factor of V2. The FWHM of the 
restoring beam of one of the component maps is displayed in the bottom 
left corner of the image. 



be a suitable dividing line between point and extended sources at 
all SNRs. 

The minimum component size, e aii , that can be measured is 
proportional to the FWHM of the synthesis ed beam and in versely 
proportional to the square root of the SNR (Fomalont 2006). Fig. [7] 



shows the median value of e ma j versus p ' , where p 



5*. for 

artifical point sources with SNRs ranging between 7 and 1000. It 
can be seen that the expected relation is obtained for p < 100. 
The relation starts breaking down at higher SNRs. This is proba- 
bly due to the limited level of accuracy with which we are able to 
simulate data, produce maps and recover source parameters. For in- 
stance, the maps have a finite pixel resolution which will affect the 
accuracy with which angular sizes can be measured. In any case, 
calibration errors are the dominant source of error for SNRs > 20. 



These results indicate that f} 



''mail' 



ought to be a useful 



indicator of source extension. Fig. [8] shows the distribution of f3 
for the point sources with SNRs of approximately 10. Fewer than 
1.2 per cent of sources have f3 > 3.0. In addition, the distributions 
of f3 for sources with SNRs ranging between 7 and 100 were inves- 
tigated. In each case, less than 1.5 per cent of sources were found 
to have ft > 3.0. As noted above, at high SNRs, calibration errors 
will become the main source of uncertainty. As a result, sources 
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Figure 6. Measured distribution of e ma j for synthetic 100-mJy point 
sources, with noise about 10 mjy. 
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Figure 7. Median value of e ma j versus p '' 2 for synthetic point sources. 
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Figure 8. Distribution of emajP 1 ' 2 /^*™ for synthetic 100-mJy point sources, 
with noise about 10 mjy. 



Figure 9. The critical component size versus SNR for b ma j = 37.5 arcsec. 



with e ma j < 25 arcsec are never considered to be extended, no mat- 
ter how high the SNR. This value corresponds approximately to the 
critical value at a SNR of 20, above which calibration errors begin 
to dominate. On the basis of our findings, a source is classified as 
extended if e mai > e cric , where 



3.0& maj p-" 2 
25.0 arcsec 



if 3.06 raaj p-'/ 2 > 25.0 arcsec, 
otherwise. 



(4) 



Clearly, there is no sharp dividing line between point and ex- 
tended sources, and no single method for making the distinction. 
However, the above procedure for estimating source sizes is sim- 
ple and robust, and is able to deal with sources with a wide range 
of SNRs. The expected percentage of point sources which are mis- 
classified as extended in this scheme is % 1 per cent. The minimum 
angular size that can be measured is estimated as a function of both 
the telescope resolution and the SNR. The median value of £> ma j for 
each survey field varies between 32.5 and 42.6 arcsec; it is gener- 
ally lower for a field lying at a higher declination. Fig.|9]shows how 



e cr it varies with SNR for b m 



37.5 arcsec. 



The method for distinguishing between point and extended 
sources can be used to establish how best to measure the total flux 
density of a source. Clearly, S m should be used to estimate the total 
flux density of an extended source. Both 5 p k and Sj n can be used to 
describe a point source. However, at low SNR, the results obtained 
using 5 i n will be significantly less accurate given the fact that the 
source is not assumed to be point-like when measuring 5 i„. In sum- 
mary, S p k should be used to describe all sources except those that 
are classified as extended, for which 5 i„ should be used. 



8 APPLICATION OF SOURCE EXTRACTION 
TECHNIQUES TO REAL DATA 

To illustrate the source extraction procedures outlined above they 
were applied to one of the survey fields, covering an area of about 
2.5 deg 2 ; the results are presented here. The noise levels in the 
core and extended field areas are typically 50 and 100 p:Jy respec- 
tively. In total, 158 sources were detected above 5cr, of which 10 
are classed as overlapping. 



CO 




S pk /mJy 

Figure 10. Comparison of integrated flux densities with peak flux densities 
for sources detected in one of our survey fields, with a line indicating equal 
flux-density values. The brightest sources have been omitted from the plot. 
Sources classified as point-like and extended are represented as filled circles 
and open squares respectively. 



The Gaussian-fitting converged successfully for all sources 
within 25 iterations. Fig. 1101 shows how S m compares with 5 p k for 
all sources, 10 of which are classified as extended. 

No over-narrow sources were detected. However, when using 
a primary beam cut-off of 0.3 to map the data several sources were 
found to have integrated flux-densities significantly smaller than 
their peak flux-densities. This is indicative of a spurious source. 
These detections were found to be separated from other sources 
(spurious or real) by a distance much smaller than a synthesised 
beam FWHM. This is implausible and resulted from the edge of a 
component map passing through a genuine source - this resulted in 
a sharp discontinuity in pixel values (Fig.ll U. 



© 2010 RAS, MNRAS 000,[T]{9] 



1 OC Survey - Data Analysis 9 




ft i J i i i i lJ 

03 01 46 44 42 40 3B 36 34 32 

RIGHT ASCENSION (J2000) 

Figure 11. A real source with two spurious detections when mapping was 
carried out using a 0.3 cut-off of the primary beam. The l-o" map noise is 
105 /i]y. Contour levels start at ±210 //Jy and increase at each level by a fac- 
tor of V2. Recovered peak source-positions are indicated with crosses. The 
FWHM of the restoring beam of one of the component maps is displayed in 
the bottom left corner of the image. 
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9 CONCLUSIONS 

In order to investigate the high-frequency-radio sky we have ob- 
served 10 survey fields, covering as 27 deg 2 to an rms noise level 
of g .2 mjy, using the AMI LA at 1 6 GHz . In an accompanying 
paper, [AMI Consortium: Davies et alj d2010h . we present the first 
results from the survey, including the deep 15.7-GHz source count. 
Here, we have concentrated on developing techniques for produc- 
ing and analysing the survey raster maps. In particular, we have: 

(1) developed systematic, automated methods for identifying 
and characterising sources in maps with varying noise levels and 
sythesised beams. 

(2) proposed a straightforward and robust method for distin- 
guishing between point and extended sources over a wide range 
of SNRs. Our method has been tested using maps including sim- 
ulated sources and noise, and has been shown to be successful in 
identifying extended emission. 

(3) applied our techniques to real sky maps and demonstrated 
that our automated techniques are useful for identifying and chrac- 
terising complex structure (see, for example, Fig.[4]l. 
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